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We propose a many-body generalization of the Z2 topological invariant for the quantum spin 
Hall insulator, which does not rely on single-particle band structures. The invariant is derived as 
a topological obstruction that distinguishes topologically distinct many-body ground states on a 
torus. It is also expressed as a Wilson- loop of the SU(2) Berry gauge field, which is quantized due 
to time-reversal symmetry. 



A topological insulator is a quantum phase of matter 
with a bulk energy gap that cannot be deformed contin- 
uously to a trivial band insulator without going through 
a quantum phase transition [l| . The most prominent ex- 
ample is the quantum Hall (QH) state 0, H, 3- One 
physical manifestation of topological orders in QH states 
is the existence of gapless edge modes, which are robust 
against any small perturbations 0, Q . 

Topological orders in QH states are associated with 
the charge degrees of freedom of electrons and accompa- 
nied by broken time-reversal symmetry (TRS). Recently, 
Kane and Mele have established a topological classifica- 
tion in time-reversal (TR) invariant noninteracting sys- 
tems associated with the spin degrees of freedom Q. 
Considering the doubled version of the Haldane model 
for the integer QH effect Q as an example, they showed 
that there exists a new type of topological insulators, 
quantum spin Hall (QSH) insulators, which are charac- 
terized by an odd number of Kramers pairs of gapless 
edge modes. In contrast, trivial insulators have an even 
number of pairs. QSH insulators in the bulk are charac- 
terized by a Z2 number 0, i, 113, EH, [II ■ 

Although the QSH insulator has been proposed in a 
noninteracting system, the even or odd parity in the num- 
ber of Kramers pairs of edge modes is robust against weak 
interactions that respect TRS 0, [H, [ill ■ This suggests 
that the Z2 topological order in the bulk is also robust 
against those perturbations. It is the purpose of this pa- 
per to construct a many-body generalization of the Z2 
invariant that detects such bulk Z2 topology. When in- 
teractions are weak, the construction we propose coin- 
cides with (and hence is a natural generalization of) the 
known Z2 invariant in noninteracting limits. While it 
is applicable even in the presence of strong interactions, 
our scheme does not exclude the possibility that the con- 
struction of the Z2 invariant might not be unique, and 
hence the corresponding strongly interacting many-body 
state can be characterized by a multiplet of Z2 invariants. 

We put a many-body system with a fixed number of 
particles on a 2D torus with N ^ x Ly unit cells, 
each of which has m sites, where N is assumed to be 
odd. We focus on the case with m = 2 which includes 
the Kane-Mele model generalized by adding interactions. 
Consider a deformation of the Hamiltonian by threading 



FIG. 1: (a) The torus in real space where magnetic fluxes 
ax,y are threaded in the a;,y-cycles, respectively, (b) The 
space of the Hamiltonian obtained by the flux-threading. 
The doublet v^{a) (u^(a)) is well-defined in the patch A(B). 



magnetic fluxes < a^^y < 27r through the two cycles of 
the torus where ax {ay) denotes the fiux threaded along 
the x{y) directions [Fig. [1] (a)]. The space of the fluxes 
forms a torus, which we call to distinguish it from the 
torus in real space [Fig. [T] (b)]. We consider a 27V- 
particle ground state. If there exists a finite energy gap 
above the ground state at a = 0, there is a gap every- 
where in for a large enough system, as the effects of 
the fluxes becomes smaller for larger system[l^. There- 
fore there exists a unique 2A^-particle ground state l^*") 
with a finite energy gap at each point a € . 

The QH states are characterized by a topological ob- 
struction of the U(l) bundle on T^, i.e., inability to 
smoothly define the ground state over T^, which is sig- 
naled by a U(l) Chern number [3, [l3|- On the other 
hand, in the QSH effect, the U(l) Chern number always 
vanishes because of TRS. However, a non-trivial sub- 
structure can still be hidden in the ground state. For ex- 
ample, when the z-component of spin (S^) is conserved, 
the Chern numbers Ch|.| for up and down spins can be 
defined separately, which are related by Ch| = — Ch| be- 
cause of TRS. The parity (even or odd) of Ch^ defines 
the Z2 invariant in this limiting situation. 

To probe such a non-trivial substructure of l^""), we 
consider the following iV-electron states obtained by cre- 
ating N holes in the ground state. 
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where cl^g/c^^g creates/annihilates an electron with 
spin s at a site t = 1,2 within a unit ceU la- 

beled by r. If d{a) :— det(l, s|l, s') ^ 0, the vector 
space spanned by the states |l,t (a)) and |1, j (a)) is 
two-dimensional, and a doublet of orthonormal states 
v{a) = (|l(a)), |2(a))) can be obtained from suitable lin- 
ear combinations, \n{a)) = ^^^-^ \^Uns{oi)\l, s{a)) with 
n = 1,2. Uns{a) is chosen so that {n{a)\m{a)) — 6n,m- 
From the doublet, we can define a projection operator 
P(a) = |l(a))(l(a)| + |2(a))(2(a)|. 

In the case of non-interacting electrons (such as the 
Kane-Mele model), the doublet v{a) can be constructed 
from the Bloch wavefunctions, and the specific way of cre- 
ating the doublet, e.g., choice of the sublattice we made 
(t — 1) at which we remove electrons, does not mat- 
ter, i.e., different choices of sublattice lead to the same 
projection operator. When we perturb the system by 
weak interactions, since d{a) cannot vanish abruptly, the 
above construction of the doublet v{a) smoothly inter- 
polates non-interacting and interacting cases. However, 
in general, d{a) can vanish at some points in where 
the two states |l,s(a)) with s =1,1 fail to span a 2D 
vector space. Since d{a) is a real function defined on the 
2D space T^, d{a) can vanish either on lines or points in 
T^. d{a) cannot vanish on a 2D submanifold of un- 
less d{a) is identically zero since the ground state l^""), 
and hence d{a) as well, are analytic on T^. If ^(0;) van- 
ishes at all points in T^, we can perturb our Hamiltonian 
slightly without closing the gap so that d{a) becomes 
nonzero everywhere except for some lines and points. 
Since the dimension of the submanifold where d{a) van- 
ishes is smaller than 2, we can analytically continue the 
2D projection operator P{a) to the whole generically. 
One may think that such a smooth continuation may not 
work because there can be points ai at which one (or a 
linear combination) of the doublet \l,s{a)) vanishes ei- 
ther due to a scalar vortex (winding in the phase of the 
coefficient of a quantum state with a vanishing norm at 
the center), i.e., |1, s{a)) = [{a^ - a^x) - a.iy)]\'4't) 

or due to a vector vortex (winding between two quan- 
tum states with a vanishing norm at the center), i.e., 
|l,s(a)) = (a^ - aix)\ipii) + {ay - aiy)|i/'j2) where 
and \ipia) are well-defined states with {ipii\^i2) = 0. For 
scalar vortices, P{a) can be smoothly extended because 
the projection operator is insensitive to the overall phase. 
On the other hand, for vector vortices P{a) cannot be 
smoothly defined. However, one can generically avoid 
the occurrence of the latter points by adding a small 
TR symmetric perturbation to the Hamiltonian. If we 
add a perturbation, the points where the norm vanishes 
will disappear because the state will generically change 
as \l,s{a)) |l,s(a)) + A{a)\ijl) where A{a) at 
a = ai and is generically independent with {ipia) 
and |1, —s{ai)). In this way, the 2D projection operator 
P{a) can be smoothly defined in the whole T^. 

We now choose v{a) such that it satisfies the TRS con- 
dition, v{—a){ia2) = &v{a), where O is the TR operator 
p^ . Note that when acting on v{a), which consists of 



TV-electron states^ 8^ = —1. While this is always possi- 
ble locally since QP{a)Q^^ ~ P{—a), there is no guar- 
antee that such basis can be defined globally. There- 
fore, we first divide into two patches {A and B) as 
shown in Fig. [T] (b) , and see whether we can merge them 
into a single patch or not. In each patch, we have a 
doublet vP{a) = (|F(a)), |2P(a))) (p = A,B), satisfy- 
ing (nP{a)\mP{a)) = Snm, vP{-a){ia2) = QvP{a), and 
Pia) = \Via)){Via)\ + \2Pia)){2P{a)\. 




FIG. 2: Two topologically different paths the SU(2) matrix 
g{a) can take as the angle ip in An B changes from to n. 
Here the "sphere" represents the space of SU(2) matrices, S^. 
The path in (a) is contractible while the path in (b) is not 
because <?(— a) should be at the opposite point of the g{a) in 

Since v^{a) and v^{a) span the same 2D Hilbert space 
a,t a e An B, they are related by v^{a) = v^{a)M{a), 
where M{a) is a U(2) matrix (transition function). Us- 
ing the TRS condition, we obtain a relation between the 
transition function at ~a and a, M{—a) = a2M*(a)a2- 
If we decompose M{a) into a U(l) phase 9{a) and an 
SU(2) matrix g{a), we have the two possibilities, 

a) 9{—a) = —0(a) and g{—a) ~ g{a), 

b) e{-a) = -61(a) + 7r and g{-a) = -g{a). (2) 

In the TR invariant system, the U(l) phase cannot have a 
non-trivial winding. However, the SU(2) matrix has two 
topologically distinctive configurations. In the case a), 
we can deform a ground state continuously such that the 
trajectory of g{a) as a function of if, where ip parameter- 
izes the overlapping region AnB which has the topology 
of S^, is contracted to a point in the space of SU(2) ma- 
trices as shown in Fig. [5] (a). On the other hand, in the 
case b), the trajectory of g{a) is not contractible because 
17(0;) = —g{—a) as is shown in Fig. [5] (b). Therefore, a 
state in the class b) cannot be continuously deformed to 
a state in class a). The relative sign between g{a) and 
g{—a) is the Z2 topological invariant. If there is no in- 
teraction, a plays the role of momenta of Bloch states 
and the topological invariant reduces to the existing Z2 
invariant [9|. In the present scheme, the Z2 invariant can 
be generalized to interacting cases. Since the general- 
ized Z2 invariant is quantized, it cannot change abruptly 
upon turning on interactions. 

Then when can the Z2 invariant change, as we tune 
some parameters (other than a) of the Hamiltonian, 
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along a path in the space of Hamiltonians? The Zi in- 
variant is well-defined if the following three conditions 
are satisfied: (1) there is TRS, (2) the ground state 
is uniquely defined, and (3) P(ol) is well-defined at all 
points in T^. Accordingly, there are three possibilities 
where the Z2 invariant may change. The two obvious 
cases are a breaking of TRS and a phase transition. The 
last possibility is the case where Pioi) becomes ill-defined 
without encountering either of the former two situations 
either because d{a) identically vanishes or because there 
occurs vector vortices. However, if P(a) is not well- 
defined at a point on the path in the space of Hamil- 
tonians, one can detour the point by slightly modifying 
the path so that Pia) becomes well-defined at all points 
along a new path. This is always possible because the 
measure of the points where identically vanishes or 
vector vortices occur is zero in the space of TR symmetric 
Hamiltonians. Namely, the codimension of the subspace 
where P(pt) becomes ill-defined is infinite. This guar- 
antees that the space of well-defined Pia) is connected 
as long as there is no phase transition. Since the quan- 
tized Zi invariant is well-defined along the new path, 
the invariant should be the same for the two end points. 
Therefore, the Zi invariant can change only when there 
is TRS breaking or a phase transition between the two 
points. Thus, two quantum states with different Zi in- 
variants are always topologically distinct. On the other 
hand, if they have the same Z2 invariant, it is generically 
unclear if they are adiabatically connected or not. From 
this point of view, it would be interesting if there is a 
many-body state which has a non-trivial Zi invariant, 
and yet is separated from the known non-interacting Z2 
insulator by a phase transition. 

The non-trivial Zi topological order amounts to our 
inability to define the doublet globally while maintain- 
ing the TRS condition v(—OL)(ia'2) = Qv(ol). If we relax 
this condition, the doublet can always be defined in a 
single patch. For a trivial case, this is true by definition, 
whereas for a non-trivial Z2 insulator, this can be done by 
defining v{a) = v^ia) for a e A and v{pi) = {a)M(oi) 
for a e S where M (a) = M {a) for a G A n B and 
M(a) = 1 at a = 0. A smooth function M(oi) which sat- 
isfies the above conditions can always be found because 
Hi(S'"^) = 0, where Hi is the first homotopy group and 
is the space of the SU(2) group. The price to pay for mak- 
ing such a single patch is that we inevitably lose the TRS 
condition relating via) and «(—«), v{—a)(ia2) = 8w(q;), 
if the Zi invariant is non-trivial. Instead, they are re- 
lated by v(^^oi)w(ct) = G)f(a), where i/j(a) is a gen- 
eral U(2) matrix. Since 0^ = —1 w{pi) must satisfy 
uF {a) = —w{—a). Especially, at the four TR symmetric 
points, (0,0), (7r,0), (7r,7r) and (0,7r), w{a) is antisym- 
metric. From the two orthonormal states |'t-(q!)) with 
n = 1, 2 at each a, we can define a U(2) Berry gauge field, 
A™"(a)da'^ = {m{a)\d\n{a)). The U(2) gauge field can 
uniquely be decomposed into the U(l) (a") and SU(2) 
(a^) parts as = a°/(2«) -|- • <t/(2z). 



The condition that sews local frames at a and —a to- 
gether naturally induces a constraint on the Berry gauge 
field configuration, 

A^{—a) = w{a)A^{a)w^ {a) — w{a)d^w\a). (3) 



Accordingly, the U(l) and SU(2) gauge fields are con- 
strained as a° (— a) = 0^(0;) — 29^^(0;); and a^(— a) • 
(J 1(21) = ■ "10(0)0^ (a) I — w(a)d ^^w^ ia) , 

where we decomposed w{ol) into the U(l) (e*^) and SU(2) 
{w) parts, w(a) = e"'(")w(a). (Note that this decompo- 
sition has a global sign ambiguity, which will not affect 
the following discussions.) At the TR symmetric points, 
w(a) is equal to iai up to sign, w(ol) ~ Pf [w(q;)] x z(T2, 
where Pf Iw] is the Pfaffian of w. 




FIG. 3: Four time-reversal invariant loops €1,2,3,4 on T^, and 
a loop d{T^/2) that encloses the half of T^. 

Following Ref. [Ill, a natural object to consider, which 
is invariant under a-dependent U(2) transformations of 
the local frame, is an SU(2) Wilson loop 



W\C] 



— TrPexp 



(4) 



where P represents the path ordering and C is a closed 
loop on the a plane. An essential observation is that 
for loops C that are invariant under TRS (i.e., loops 
that are mapped onto themselves by TRS up to orien- 
tation), the sewing condition ^ quantizes the SU(2) 
Wilson loops, W[C] = ±1. For example, for a straight 
loop Ci running from (-7r,0) to (7r,0) (Fig. [3]), W[Ci] 
Pf [tI;(7r,0)]Pf [w(0,0)]. 

Although being invariant under a gauge transforma- 
tion induced by a unitary transformation of the local 
frame at a, each Wilson loop is not invariant under re- 
definition of unit cells. For example, when there are only 
two orbitals (t = 1,2) in each unit cell, a spin-dependent 
redefinition of unit cell Cr,i,s — > Cr,i,s, Cr,i,s Cr+spi,2,s, 
p G Z flips the sign of M^pi] and M^[C2] when p is odd, 
where x is a unit lattice translation vector in x direction. 
In other word, this transformation inserts a half unit of 
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the SU(2) flux along the ay-direction. When there is 
translation invariance, the choice of unit cells is arbi- 
trary, whereas if we introduce a boundary, the choice of 
unit cell should be consistent with the location of the 
boundary. In this sense, the (quantized) value of each 
Wilson loop matters when we terminate the system. 

We can construct, however, from two parallel Wilson 
loops (e.g., Ci and C2), an invariant which is left un- 
changed by the half-flux insertion in ax and ay direc- 
tions, since the effect of the half- flux insertion cancels. 
We are thus led to consider a Wilson loop that runs 
counterclockwise around the boundary of the half of Tj, 
T2/2 := (-TT,^] X [0,7r], say (Fig.©: 

W[din/2)]:={-l)^= n PfKfc)]- (5) 

fe=(0,0),(7r,7r), 
(7r,0),(0,7r) 

The Z2 number A distinguishes trivial (A = 0) and non- 
trivial (A = 1) insulators. The Kane and Mele modelQ 
is an explicit example of the latter case. 

The connection between the gluing and the gauge pic- 
tures can be established for the non-interacting case with 
conserved . In this case, the ground state can be writ- 
ten as l^f") = 1^1") X Iv^J) and the doublet becomes 
v{a) = (l^fp, l^*")). The Z2 invariant obtained in the 
gluing picture is nothing but (— 1)'-'^t = (_i)Ch^_ On the 



other hand, the quantized SU(2) Wilson loop also gives 
W[d{T^/2)] = (-l)ChT. The Z2 invariant is independent 
of the choice of sublattice in Eq. H]) because annihilating 
electrons of spin s results in the same state of l^"" ,) ir- 
respective of the choice of sublattice. Since the quantized 
Z2 number cannot change abruptly as we turn on inter- 
actions or disorders, it is the only topological invariant in 
weakly interacting TRS systems. Thus the Z2 invariant 
is a natural generalization of the parity (even/odd) of the 
Chern number, to the case without the conservation 
in which case the Chern number cannot be defined. 

Finally, a few remarks are in order. First, if we add 
some perturbations which enlarge the size of the unit 
cell from to = 2 to m > 2, we can still remove half of 
the total electrons at the original location to calculate 
the Z2 invariant. This allows one to characterize the 
topological order of systems which have both interactions 
and disorder Second, the formahsm developed so far 
can be extended to three dimensions similarly [l^, [l^, 

m. 
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